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Quasi-particle picture in a magnetic field is pursued for dynamical spin and charge correlation 
functions of the one-dimensional supersymmetric t-J model with inverse- square interaction. With 
use of exact diagonalization and the asymptotic Bethe-ansatz equations for finite systems, ex- 
citation contents of relevant excited states are identified which are valid in the thermodynamic 
limit. The excitation contents are composed of spinous, antispinons, holons and antiholons 
obeying fractional statistics. Both longitudinal and transverse components of the dynamical 
spin structure factor are independent of the electron density in the region where only quasi- 
particles with spin degrees of freedom (spinous and antispinons) contribute. The dynamical 
charge structure factor does not depend on the spin-polarization density in the region where 
only quasi-particles with charge (holons and antiholons) are excited. These features indicate the 
strong spin-charge separation in dynamics, refiecting the high symmetry of the model. 

KEYWORDS: dynamical spin structure factor, dynamical charge structure factor, supersymmetric t-J model, 
inverse-square interaction, magnetic field, spin-charge separation, exact diagonalization, recursion 
method 



§1. Introduction 

The concept of spin-charge separation is no longer 
a purely academic idea. In various low-dimensional 
electron systems, angle-resolved photoemission measure- 
ments have revealed the existence of spin-charge separa- 
tion. In quaRirOne-dimensionaLQD) Mott insulators such 
as SrCu02a'B' and Sr2Cu03,BQ' angle-resolved photoe- 
mission data show two distinct dispersions in the energy- 
momentum plane. The two correspond to the spinon 
band and the holon band, respectively. This is an indica- 
tion of the spin-charge separation. It has been reported 
that the spin-charge SjSparation in Sr2Cu03 appears also 
in dielectric response Ja-* 

Recently magnetic-field-induced dimensional crossovjer 
has been observed in the normal state of YBa2Cu408l3-' 
The CuO chains of this material play the critical role 
in magnetoconductivity. The following question may 
arise: When a ID electron system is under a magnetic 
field, how is the spin-charge separation affected by the 
field? The magnetic field controls spin degrees of free- 
dom, while the hole doping controls charge degrees of 
freedom. However, realistic situation should be rather 
complex due to convolution of both degrees of freedom 
in a strongly correlated system. 

It is well-known that interacting electron systems in 
one dimension are described by the Tomonaga-Luttinger 
liquid,LlEP but the information is restricted to that in the 
long- wavelength and low-energy limit. One can adopt 
the Hubbard model and the supersymmetric t-J model 
as Bethe-ansatz solvable models. With use of the exact 
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wave function of the C/ — > cxd Hubbard model, critical ex- 
ponents and global features of static spin structi. 
tors have been investigated in a magnetic field.l 
The exponents depend only on the magnetization per 
electron and are independent of the electron density. For 
the supersymmetric t-J model with nearest-neighbor in- 
teraction, the exponents depend both on the magnetiza- 
tion and the electron density.H In order to obtain deeper 
understanding of spin-charge coupling, it is imperative to 
derive theoretically not only the critical exponents but 
also the overall spectral weight in dynamics. 

The supersymmetric f- J model with inverse-square in- 
teractionllj) reveals the spin-charge separation in a typi- 
cal form. The spin velocity is independent of the electron 
density n, and the charge velocity is independent of the 
spin-polarization density ifix3 At low temperature, the 
spin susceptibility is independent af p, and the charge 
susceptibility is independent of m.t§E3' It has been found 
from the exact thermodynamics that elementary excita- 
tions consist of four types of quasiLparticles: spinous, 
antispinons, holons and antiholons .ES' In zero magnetic 
field, the regions of nonvanishing spectral weight in the 
momentum-energy plane, which are called the "compact 
support" , have been obtained for various dynamical cor- 
relation functions at zero temperature.EZP The excitation 
contents are composed of spinous, holons and antiholons. 
In particular, the weight of the dynamical spin structure 
factor does not depend on n in the region where only two 



PC) 



ling with n = 1, i.e., in 
magnetic-field dynam- 



spinons contribute.EEP At haJf- 
the Haldane-Shastry modeljlj' 

ics has been investigated. Ell' So far, study of full dynam- 
ics are lacking at arbitrary filling and spin-polarization 
density. 
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In this paper, we extend our previous workll3^ to the 
case of finite magnetic field for the l/r^-type supersym- 
metric t-J model. We investigate how the spin-charge 
separation appears in the spin and charge dynamics in 
nonzero field. To this end, exact diagonalization up to 
16 sites is used. In order to analyze the results of dynam- 
ical quantities, we employ the asymptotic Rpthp-anpaty 
equationsot^ and the skew Young diagram.tSBEJ' In 
addition, we examine the static structure factors in the 
1/r^ model as well as the i- J model with realistic param- 
eters. 

This paper is organized as follows: Section 2 describes 
the model and definitions of dynamical quantities. In 
§3 we derive the supports for dynamical spin and charge 
structure factors. In §4 we discuss the spectral weights 
themselves, and clarify manifestation of the strong spin- 
charge separation in dynamics. In §5 we show the results 
of the static structure factors. Finally we give the sum- 
mary of this study in §6. A brief account of this work 
has been presented in rcf. 25. 

§2. Model and Dynamical Structure Factors 

We consider the following Hamiltonian: 



n 



N 

E 



+JU 



T.^h 



CirrCja -\- 



Cj^Ci, 



■ Sj - -^n^"^: 



N 

1=1 



(2.1) 



with even number N of sites and average electron number 
n per site. Here Ci^ = Ci„(l — ni^^„) is the annihilation 
operator of an electron with spin a at site i with the 
constraint of no double occupation. The Zeeman term 
couples the system to a uniform magnetic field H in the z 
direction. A dimensionless parameter h is defined as h = 
gfinH/t. We assume the periodic boundary condition. 



Then lij is given by /i- 



[iN/Tr)smiTr{i~j)/N)] 



for the 1/r type, and by lij — Sj^i+i for the nearest- 
neighbor type. In this paper we mainly discuss the 1/r^- 
type supersymmetric i- J model {J/t = 2). 

The dynamical spin and charge structure factors are 
given by 



N{q,uj) 



E 
E 



(i.|s^|0)| 6iLO~E,+Eo), (2.2) 
{i^\ng\0)\^S{uj-E^ + Eo), (2.3) 



where s^ = A^^-^/^ ^^ s"e~"?^ (a = x,y,z), Uq = 
^-1/2 J2g{n£ — f7,)e~'*^, and \iy) denotes an eigenstate of 
Ti. with energy E^ {Eq being the ground-state energy). 
The transverse component S^^{q,uj) = S'^'^lq^tu) is de- 
composed into two parts: 



5^^(q,w) = 5^^(q,tj) 



> 



-(g,c.) + 5^ 



(2.4) 

We calculate S ~^{q,Lj) and S~^ {Qt^) instead of 
S^^{q, ui) and ^^^(q, oj). These dynamical quantities can 
be written in the form of a continued fraction.E^' The 
coefficients of the continued fraction are obtained from 



the Lanczos algorithm. This method is called the re- 
cursion method. We truncate the continued fraction be- 
tween 100 and 200 iterations and take the Lorentzian 
width of C'(10~''^i). In our calculation, energy levels and 
intensities have accuracies of 7 digits and 4-5 digits, re- 
spectively. 

We assume m > without loss of generality. In the 
thermodynamic limit, fh is determined by h as follows: 



^=( l-v/T^W^, f"^;-^-"- (2.5) 

I n tor h> he 

where he = n(2 — njn^ /2S3 For finite systems, ifi — 
{N -Q - 2M)/N and n = [N - Q)/N. Here Q and 
M denote the numbers of holes and down-spins, respec- 
tively. For n < 1 we consider the case of Q = even and 
M ~ odd under the periodic boundary condition. 

§3. Spectra of Quasi-Particles and Supports 

In this section we determine the regions of nonvan- 
ishing spectral weight in the energy-momentum plane 
(i.e., the compact supports) for various dynamical cor- 
relation functions. To achieve this goal, the following 
points should be clarified: 

• Dispersion relations of the relevant quasi-particles. 

• Excitation contents of the excited states. 

• Selection rules for the excited states. 
Let us begin with the dispersion relations. 

3. 1 Dispersion relations of quasi-particles 

The spectra of quasi-particles in the model are most 
easily obtai|a£d_lm.the asymptotic Bethe ansatz method. 
It is knownEjLS'Ej) that the method gives a part of the 
exact eigenstates corresponding to the Yangian highest- 
weight states (YHWS), but that other eigenstates are 
missed. Since the YHWS contain all kinds of quasi- 
particles, the asymptotic Bethe ansatz is sufficient to 
derive the dispersion relations. The Yangian multiplet 
structure of eigenstates is then derived by the skew 
Young diagram. H 

We write the energy of a spinon with spin a as e^, while 
that of an antispinon as £§. Note that an antispinon 
can have only a down spin with m > 0. The energies 
of a holon and an antiholon are written as eh and ejj, 
respectively. Details of the derivation of the spectra are 
given in Appendix. The final results are given by 



el/t = -k{k-TT)--, 


(3.1) 


ei/t = -k{k-n)+'^, 


(3.2) 


4/t = ^Hk - 27r) + h, 


(3.3) 


^'^/^=('-2) +12' 


(3.4) 


1 ., ^2 7r2 


(3.5) 



where h — m(2 — m)7r^/2. 

We note that only a part of the momentum space is 
relevant as shown in Fig. 1. For classification of ex- 
cited states, it is convenient to define the right and left 
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branches for each kind of quasi-particles. For example, 
the right [left] up and down spinons are allowed only 
in the range ffnT/2 < k < nTr/2 [(1 — n/2)n < k < 
(1 — fh/2)7r]; the right [left] holons are allowed only for 
(1 - n/2)n < fc < (1 - fh/2)Tr [fhn/2 < k < nTr/2]. 
The range of momentum for the right [left] antispinons 
is (2 — fh)TT < k < 27r [0 < k < rhir], while the antiholons 
propagate in the region fn: < k < {2 — n)n . 

By shifting the origins of momentum and energy for 
each kind of quasi-particles, we can write the dispersion 
relations for right (left) spinons esnjLji right (left) anti- 
spinons Csn/L), right (left) holons EhR/L, and antiholons ejj 
in the following form: 



^Sr(l) 



/t = -k{k T «°) 



1 



1 r 



ehA=o ( 



(3.6) 
(3.7) 
(3.8) 
(3.9) 
n)Tr. In the shifted 



where v^ — (1 ~ fa)i: and v^ — {\ 
momentum space, the right [left] spinons and holons are 
allowed in the range < fc < {n — 7fi)TT/2 [—{n — ifi)TT/2 < 
k < 0]. The allowed region of momentum for the right 
(left) antispinons is < fc < ffn: {—frn: < fc < 0), while 
the antiholons propagate in the region — w° < fc < w°. 
In the case of rfi = 0, these dispersion relations and the 
range of jnomentum reduce to those given by Ha and 
Haldane.0) 

3. 2 Determination of the supports 

The remaining task to determine the supports is to 
derive excitation contents and selection rules for each 
dynamical correlation function. 

First, let us consider excitation contents. For a 
given magnetic field, the magnetization of the ground 
state is fixed. In the subspace of relevant excited 
states, we use the solution of total energy and to- 
tal mpmentum from the asymptotic Bethe-ansatz equa- 
tions.oE3' We then obtain the motif corresponding 
to the ground state and excited states with finite 
intensity obtained by the recursion method. The 
motif of the ground state with finite magnetization 
is given by 00.. .00 | 1010.. .10 | 11. ..11 1 01. ..0101 1 00.. .00. 

I ,11 III IV . . Y 

This ground-state motif turns out to be divided into 

five segments (TV). The segment III can be regarded 
as a holon condensate, and the segments I and V can 
be regarded as a up-spinon condensate. As a pertur- 
bation from the magnetized ground state, spinons (00) 
and holons (11) are excited in the segments II and/or 
IV; Antiholons (0) are created in the segment III, and 
antispinons (1) are created in the segments I and/or V. 
Spinons and holons have semionic statistics, while anti- 
spinons and antiholons have Bose statistics. 

Let S^ and C denote the z-components of spin and the 
charge of each quasi-particle, respectively. It is natural 
to assign {S^,C) = (1/2, 0) to the up spinon, (-1/2, 0) to 
the down spinon, (—1, 0) to the down antispipoB_(0, +e) 
to the holon and (0, -2e) to the antiholon.lia'B'k§) How- 



ever, this assignment leads to cases where some states 
with excitation contents decided above do not fulfill the 
spin and/or charge conservations for the excited states. 
Moreover, those states do not fulfill the statistics of 
quasi-particles, either .J^his situation has occurred also 
in zero magnetic field. Ej-' 

In ref. 18, we introduced the skew Young diagram in 
the l/r^-type supersymmetric t-J model. By using this 
diagram, one can reproduce the supermultiplet structure, 
i.e., nontrivial degeneracies in the excitation spectrum. 
In addition, excitation contents for a finite system can be 
extracted by the skew Young diagram. This idea can be 
expanded into the magnetized case. We first determine 
the five segments of a diagram from the corresponding 
motif. Note that the segment HI consists of Q boxes, 
and that each of the segments I and V consists of {N — 
Q — 2M)/2 boxes. Next, we put an index for the internal 
degrees of freedom on each box of a skew Young diagram, 
according to the rules presented in ref. 18. We represent 
an up-spin by 1, a down-spin by 2, and a hole by o. One 
can read quasi-particles from a diagram of an excited 
state as follows: 

(i) We regard o in the segments II and IV as a holon. 

(ii) We regard a row of a 1-2 pair in the segment III 
as an antiholon. 

(iii) In the segments II and IV, single 1 or 2 in a col- 
umn, both of which do not make a pair, is regarded as a 
spinon; it does not matter whether the column include o 
or not. Single 1 corresponds to an up-spinon, and single 
2 to a down spinon. 

(iv) A column of a 1-2 pair in the segments I and V is 
regarded as an antispinon. 

Because an antiholon or an antispinon consists of two 
boxes, there exist cases where only a part of a 1-2 pair is 
present in the segment defined by (ii) and (iv). We shall 
call such a part of the 1-2 pair a "half-antiholon (h*)" for 
the former, and a "half-antispinon (swr-i)" for the lat- 
ter. The half-antiholon [half-antispinon] is assumed to 
have {S'^^C) = (0, — e) [(—1/2,0)] and semionic statis- 
tics. Thereby the excitation contents from the diagram 
satisfies both charge and spin conservations and the sta- 
tistical feature. 

Let us give a example of the skew Young dia- 
gram. For S^ [It^) from the initial state with 
{N,Q,M) — (12,2,3), the excitation energy with 
[ql 71,10 It) — (0.5,2.05617) has finite intensity. Its motif 
is 010|10|110|10|100, and the corresponding skew Young 
diagram is shown in Fig. 2. The excitation content reads 
Sl -|- (sL,hL) -I- ii* + s'^- Then the change of the z- 
component of spin becomes —1 and there is no change 
of charge, which is consistent with the changes between 
the initial state and the final one. 

As the number of sites increases with electron density 
and spin-polarization density kept constant, the effect 
of the edges in the holon condensate and that in the 
up-spinon condensate should be negligible. This means 
that both the half-antiholon and the half-antispinon have 
zero energy. In fact, for N{q,uj) at arbitrary filling and 
zero field, the energy levels with contents h* + hR are 
almost along the holon dispersion for < g < n7r/2 
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in the compact support JiS-* For S^^{q,uj) at n = 1 and 
arbitrary magnetization, the energy levels with contents 
s + sj!;^ are almost along the spinon dispersion forJl< q < 
(1 — m)7r in the support by Talstra and HaldaneO Thus 
half-antiholons and half-antispinons should not survive 
in the thermodynamic limit. In this way, the excitation 
contents in the thermodynamic limit are identified as 
shown in Table I. The z-component of spin (1/2 or —1/2) 
of spinous is assigned so that the spin conservation is 
satisfied. 

When we draw the support using the excitation con- 
tents (quasi-particles are identified as free particles), the 
estimated region is much larger than the region where 
the poles with finite intensity are present. This implies 
the importance of selection rules. Ha and HaldanetlP 
found the empirical rule in zero magnetic field: for a 
given spinon- holon pair (sr, h^), |usj,| > \vhj^\, and the 
same to the left-going pair. Here v is given by de/dq. 
We find that this rule applies also to the case of nonzero 
magnetic field. Then the estimated region is found to 
be the same as the region expected from the recursion 
method. 

3.3 Characteristic features of the supports 

Figures 3-9 show the regions of compact support for 
S-+{q,Uj), 5^^((7,w), S+-{q,uj) and Niq^Lu) with finite 
magnetization. The solid lines mean either the spectrum 
of a quasi-particle or that of a pair moving with the same 
velocity. The presence of other quasi-particles excited at 
an edge of each band causes several lines to appear with 
particular shifts in momentum and energy. Outermost 
dispersion lines correspond to the boundary of the com- 
pact support. Namely, there is no intensity outside the 
outermost lines. 

The number of zero modes is increased by the external 
magnetic field. Namely for S ^{q,uj) and S'"' {q,uj), 
w = is allowed at g = tott, hn, (2 — h)'K and (2 — m)7r. 
For S^^{q,u!), on the other hand, w = is allowed at 
q — [27r], 2m7r, {h — m)7r [— w" — v°], (2 — n — m)7r 
[— v'^ + v^], {n + ?7i)7r, (2 — n + ?7i)7r and (2 — 2?7i)7r. 
For N{q,Lj), w = is allowed at g = [27r], (2 - 2n)7r, 
(n — m)7r, (2 — n — r?i)7r, (n + 171)7:, (2 — n-l-77i)7r and 2n7r. 

For S ^{q,u!), the excitation contents, is the same as 
those for S{q,iLj) in zero magnetic fieldtZl' (see Table I). 
Note that no antispinons contribute to S ^{q,u!). At 
arbitrary filling, the support turns out to be essentially 
a squeezed version of S{q, uj) in zero field. 

For S^^{q,Lj) and N{q,Lj), the region of (sl, hL) + h 
+ (liRj Sr) shrinks and splits into three directions by a 
magnetic field: the two of them is the parallel shift by 
zbmTT along the momentum axis (horizontal axis), and 
the remaining one is the parallel shift by +h along the 
energy axis (vertical axis). This is due to the presence 
of the antispinon. 

Finally we note that the support of S^ (q,^) over 
h agrees with that of S^^{q,Lu) for the same filling and 
magnetization. 

§4. Spectral Weights 

In this section we focus on spectral weights in the com- 
pact supports. 



4.1 S-+{q,co) 

First of all, we discuss the results of S ^{q,uj) at 
n — 1. Figures 3(a) and 3(b) show the results in the 16- 
site chain for m = 0.25 (10 up-spins and 6 down-spins) 
and 771 — 0.5 (12 up-spins and 4 down-spins), respec- 
tively. With use oi h = fh{2 — m)7r^/2, the values of the 
magnetic field are taken as h = 2.15898 for ffi = 0.25 
and h = 3.70110 for m — 0.5. We find that for finite sys- 
tems the momentum g, the excitation energy uji, and the 
form factor MJ have strong correspondence with those 
for S{q,u!) in zero field. Namely, for < C2 <-ci-4C M, 
with ci and C2 being integers, it is known thatc^^' 



27r , , 

q^rmr+ — {ci + C2) , 



""^^iCi^ [^M{c,+C2)-2{cl + cl) 



M?, 



c\ -~ C2J , 
2ci - 2c2 



(4.1) 



(4.2) 



2ci - 

ci-l 

^ n 



1 



1 1 

~ ' 2M - 2c2 - 1 
2i 2M - 2i 



2i-l 2M-2i-l^ 

for C2 < ci - 1, 



(4.3) 



1 



(2ci - 1) (2M 



2ci + 1) 
for C2 



ci-1. 



Here M denotes the number of down-spins in the ini- 
tial state. The number of poles with finite intensity is 
M{M + l)/2. We have checked the validity of eqs. (4.1)- 
(4.3) by comparison with numerical results up to A^ = 16. 
Taking the thermodynamic limit, we obtain the ex- 
plicit formula 



S-+{q,Lu) 

_ 1 e(eu(g)-co)e(^-eL-(g))9(c^-£L+(g)) 
2 ^(w-eL-(q))(w-eL+(g)) 

where Q{lu) is the step function, and 

eL-(g) = t{q - ffnT){TT - q), 

<^L+{q) = t{q -271 + m7r)(7r - q), 



(4.4) 



(4.5) 
(4.6) 



<^u(g) == x(9 - 27r + m7r)(m7r - g). (4.7) 

In the zero- field limit (m = 0), this expression reduces 
to the result obtained by Haldane and Zirnbauer.EiP 

Next we consider the case of n < 1. Figures 4(a) and 
4(b) show the results of S ''(g, w) in the 16-site chain 
with n = 0.875 (2 holes) for to = 0.25 (9 up-spins and 
5 down-spins) and to = 0.5 (11 up-spins and 3 down- 
spins), respectively. Like the zero-field case (to = 0)}3 
at fixed to(> 0) energy levels and intensities in the two- 
spinon region (i.e., the hatched area in Fig. 4) agree with 
those for n — \ within the numerical accuracy. We have 
checked this fact for iV < 16 with various n and to. 
From this fact, it is conjectured that analytical expres- 
sion of S ^(g, w) in the 2sr and 2sl is identical with eq. 
(4.4). This feature is an indication of the strong spin- 
charge separation, which should be due to supersym- 
metric Yangian symmetry of the present model. For the 
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nearest-neighbor supersymmetric <- J model, such strong 
separation does not occur. This statement holds also for 
S^^{q,uj), S^ (q,^) and N{q,cu) to be discussed later. 

For finite systems, the momentum, the excitation en- 
ergy and the form factor are expressed as follows: 

q^rmr + ^^c, (4.8) 



N 



OJp 



MQ,M), 



M«=g,(Q,M). 



(4.9) 
(4.10) 



Here M and Q denote the numbers of down-spins and 
holes in the initial state, respectively, and c is an integer 
which satisfies < c < Q -I- 2M. We note that fi,{Q, M) 
and g^{Q, M) are independent of the system size N . This 
implies that 5* '"(g, lj) for the present model is equivalent 
to the hole propagator, where one boson is removed, for 
the SU(1,1) CalnKPro-Sutherland model with Q fermions 
and M bosons.BB) 

12 S'^^(g,w) 

Figures 5(a) and 5(b) show the results of S^^{q, oj) for 
the 16-site Haldane-Shastry model with magnetization 
m = 0.25 and 0.5, respectively. The intensity at {q,uj) — 
(0,0) is given by Nifi^/A. It is found that dominant 
intensity lies in the region surrounded by four lines: A, B, 
C and D in Fig. 5. Analytical expression of S'^^{q,Lo) in 
the region q < mn has been computed with use of strong 
coupling limit of the spin Calogero-Sutherland model. Lj' 

In the 16-site chain with n = 0.875, the results of 
S^^{q,u!) for m = 0.25 and 0.5 are shown in Figs. 6(a) 
and 6(b), respectively. The compact support reveals the 
region where only quasi-particles with spin degrees of 
freedom (i.e., spinon and antispinons) contribute. In this 
region, energy levels and intensities agree with those for 
n = 1 with the same ffi. This is an indication of the 
strong spin-charge separation in the longitudinal compo- 
nent of spin dynamics. 

4.3 S+-{q,Lu) 

In Figs. 7(a) and 7(b) we show the results of S^ {q, uj) 
in the 16-site chain at n = 1 for ?fi = 0.25 and 0.5, re- 
spectively. At g = there is a single pole; its energy level 
is Lu = h and its intensity 1(z^|5'^q|0)P = ifi. The lower 
edge over < q < ffnr is not only the boundary of the 
continuum but also the disppjsion line of an antispinon 
(i.e., magnon) excited alone.E2|) 

At a fixed n = 0.875 we show the results for the 16- 
site chain with rh = 0.25 and 0.5 in Figs. 8(a) and 8(b), 
respectively. Let us compare the present case with the 
half-filled case in the same magnetization. Although the 
structure of intensity is partially broken down by exis- 
tence of holes, the (q, a;)-region lies where the structure 
remains intact. Namely, energy levels and intensities in 
the region agree with those for fi = 1 and the same 
magnetization. This region corresponds to a case where 
holons and antiholons are never excited, in terms of the 
support in the thermodynamic limit. This fact indicates 
the strong spin-charge separation in S~^ (q,"^)- 



U N{q,Lo) 

We now turn to the charge dynamics. In rcf. 18 we 
have discussed the features of N(q, uj) in zero magnetic 
field. Figures 9(a) and 9(b) show our results for electron 
density n — 0.875 (16 sites and 2 holes) with two values 
of magnetization ffi. The following remarkable feature is 
found: in the pure h -|- 2hR region, the poles and intensi- 
ties of N{q, uj) are independent of m within the numerical 
accuracy. For q < {n — m)7r/2, in particular, this state- 
ment has been proved analytically.^ The independence 
is a manifestation of the strong spin-charge separation in 
charge dynamics. 

§5. Static Structure Factors 

Integration of the dynamical structure factor over lo 
yields the static structure factor. Figures 10-12 show 
numerical results of S^^{q), S^^{q) and N{q) for the 
inverse-square type. 

Particularly for S^'^{q), the analytic expression can be 
obtained. At n = 1, by integrating eq. (4.4) over ui, we 
obtain 

0, for < q < fhn, 

S-+{q)^{ 1, l-q/n 



2 1—771 



for 7777r < q < n. 



(5.1) 



(5.2) 



The Fourier transformation leads to 

{Sq S+) = ^:^Si ((1 - m)Trr) , 

where r is an integer which satisfies r > 1, and Si(a::) is 
the sine integral. 

For 77 < 1, numerical calculations show that S ^{q) 
are equivalent to that for n = 1 if < g < nn; S ^{q) 
for 777r < g < TT are equal to the value for q — nw. 
On the basis of these facts, we conjecture the following 
expression: 

0, for < g < 7777r, 

S- + {q) = { -2 1°S -Y~^' ^°'' m^<1< -n-^, 

1 , l~n ^ _ 
"t: iog —^ lor rnr < q < TT. 

2 1 — 777 

(5.3) 
The Fourier transformation leads to 

{S,S+) = i^ [Si ((1 - 77-7)^r) - Si ((1 - 77)^r)] . 

(5.4) 
When-ja = 0, this result reduces to the previous re- 
suIlBN) 

The relation between S ^{q) and 5*^ (g) is as follows: 



S+-{q) = m + S-+iq), 



(5.5) 



irrespective of the electron density. Therefore, using the 
relation S"=^(q) = [S^+iq) + S'+"(q)]/4, we get analytic 
expression of S^^{q): 



S--{q)^^ffl + ^S-+{q), 



(5.6) 



where 5* ^{q) is given by eq. (|5.lD or eq. (5.3). 

The results of static structure factors for the inverse- 
square type are summarized as follows: 

(i) At fixed 777, S^^{q) is independent of 77 for < q < 



Yasuhiro SAIGA and Yoshio KuRAMOTO 



[n — fri)TT. 

(ii) At fixed to, S^^{q) is independent of n for < q < 
nn. 

(iii) At fixed n, N{q) is independent of to for < q < 
{n — fnjn. 

These features are regarded as a manifestation of the 
strong spin-charge separation in static structure factors. 

What about the case of the nearest-neighbor interac- 
tion? Figures 13, 14 and 15 show the results of S^^{q), 
S^^{q) and N{q), respectively, for the nearest-neighbor 
type with two representative values of J/t. We look 
at the difference of 5'""(g) [a = z,x] between the two 
values of filling (n = 1 and 0.875), and that of N{q) 
between the two values of magnetization [fh = and 
0.25). In any case, within the special momentum range 
described in (i)-(iii), the difference is less than 6% for 
S'^^(g), 5% for S'^^(g) (g = nn not included), and 1% for 
N{q). This indicates that the ID i- J model with realistic 
parameters shares the properties of static structure fac- 
tors in the 1/r^-type supersymmetric i- J model. The in- 
dependence of static spin (charge) structure factors with 
varying filling (magnetization) spreads over considerably 
wide range of momentum, if the system is near half-filling 
and under a weak magnetic field. This behavior might 
be detected experimentally in a quasi-lD conductor. 

§6. Summary and Discussion 

We have investigated the dynamical and static prop- 
erties of the ID supersymmetric t-J model with l/r^ 
interaction in a magnetic field. We have determined 
the {q, w)-region of nonvanishing spectral weight for 
S"'°'{q,Lu) [a = x,y,z] and N(q,uj) in nonzero magnetic 
field. We have found the following features about the 
spectral weights themselves: 

(1) At fixed to, S'""((7, cj) is independent of n in the 
region where only quasi-particlcs with spin degrees of 
freedom (spinous and antispinons) contribute. 

(2) At fixed n, N{q,Lu) is independent of to in the 
region where only quasi-particles with charge (holons and 
antiholons) contribute. 

These features constitute further manifestation of the 
strong spin-charge separation in the l/r^-type supersym- 
metric t-J model. 

When a system deviates from the 1/r^-type model, ele- 
mentary excitations in a magnetic field do not corx 
to spinous and holons but to mixtures of them.Ej'tJ 
This is related to the fact that the system loses the Z2 
symmetry between the statistical properties of up- and 
dowriiSpin electrons in the presence of finite polariza- 
tioncP In fact we have checked for the polarized nearest- 
neighbor t-J model that the long-wavelength limit of 
the spin and charge excitations have a common veloc- 
ity. Thus the spin-charge separation is not realized any 
longer. On the other hand, it has been shown that 
spinous and holons together with their antiparticles span 
the complete set in the Fock space of hard-jOpp. lattice 
fermions with any filling and polarization.talla) The 
completeness is independent of whether or not spinous 
and holons form eigenstates of the model. Thus the ab- 
sence of spin-charge separation in the polarized nearest- 
neighbor t-J model means that there are residual inter- 



actions between spinous and holons. 

In contrast, the l/r^-type model maintains the Z2 in- 
variance even with finite polarization. This can be un- 
derstood from the statistical matrix given in ref. 16. The 
parallel dispersion relations of up and down spinous as 
shown in Fig. 1(a) are also consistent with the Z2 in- 
variance. As we have shown in this paper, spinous and 
holons remain elementary excitations, and are essential 
to understanding dynamics for any polarization in the 
1/r^-type supersymmetric i- J model. This property en- 
titles the model to be called a "canonical" or a kind of 
"free" Hamiltonian, which is nothing but another repre- 
sentation of the higher symmetry in the model. 
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Appendix: Derivation of Dispersion Relations 
of Quasi-Particles 

In this appendix, we derive dispersion relations of four 
types of quasi-particles (i.e., spinous, antispinons, holons 
and antiholons) for the 1/r^-type supersymmetric t-J 
model in a magnetic field. The derivation is based on 
the asymptotic Bethe-ansatz equations. 

Let us take the completely up-polarized state as the 
reference state and consider M -\- Q pseudo-particles 
which represent M electrons with spin down and Q holes. 
The asymptotic Bethe ansatz leads to the following equa- 
tionS® 



^ M+Q 



t 



p=l 



elPp) ■ 



-Q I 



1 

iV2 



N-Q 



M 



(A-1) 

where e(p) = p{p — 2tt)/2 for < p < 2tt. For later 
convenience the origin of pseudomomenta is shifted by 
TT. Then e(p) = e(p + n) = {p^ — 7r^)/2 for \p\ < n. 
When we consider the value of the physical momentum, 
we r estore the origin of pseudomomenta. Note that eq. 
(A-1) includes the Zeeman term. The pseudomomenta 
p^ are determined by the following equations: 



27r J^ = Pf_,N + TT ^ sgn(p^ - qi) 



M+Q 



-TT ^ sgn(p^ -p,,) = z{pi,)N, (A-2) 

M+Q 

27r/i = TT ^ sgn{q, -p^) = w{qi)N, (A-3) 

where J^ G [-(TV ~ M - l)/2, {N - M - l)/2] for /i = 
1, 2, • • ■ , M+Q, and /, G [-(Af-f Q)/2, {M+Q)/2] for t = 
1, 2, • • ■ , Q. The quantum numbers J^ and Ii are integers 
or half- integers, and are arranged in the ascending order. 
Before considering the dispersion relations, we must 
derive the ground-state distribution functions of pseu- 
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domomenta. For the ground state, both J^ and U are 
distributed densely and symmetrically with respect to 
zero. Namely, in the subspace with {N,Q,M), J^ = 
-{M + Q - l)/2 + ^ - 1 for /^ = 1, 2, • ■ • , M + 0, and 
h = -((5-l)/2 + ?;-lfori== 1,2,---,Q. In the large- A^ 
limit, using 2T:a{p) = dz{p)/dp and 27rp{q) =^ dw{q)/dq, 
we obtain 



ao{p) = e{Bo 



1 

2^ 



D„ 



dqpo{q)S{p - q) 

(A-4) 



dp'ao{p')S{p-p') 



po{q) = eiDo - \q\) / dpao{p)S{q - p), (A-5) 

J-Bo 

where &{p) is the step function. Here we have attached 
the subscript to quantities corresponding to the ground 
state. The cut-off parameters Bq and Dq {Bq > Dq) are 
determined by 



Thereby we get the explicit expression of ai (p) . 

(iv) We obtain the excitation energy and the momen- 
tum transfer by using ai (p) , and then derive the disper- 
sion relation. 

A.l Spinous 

In a magnetic field, the dispersion relation of up 
spinous should be different from that of down spinous 
due to the Zeeman splitting. First, in order to obtain 
the up-spinon dispersion, we consider the spin excitation 
in the subspace with (A^, Q+1,M — 1), where a down spin 
is removed. We take the numbers li {i = 1,2,- ■■ ,Q + 1) 
as 



/,; 



-(g-i)/2 + z-i. 



(A- 14) 



which means that the charge excitation (holon excita- 
tion) takes the minimum energy that vanishes in the 
thermodynamic limit. On the other hand we choose the 
configuration of J^ {p, — 1,2,- ■ ■ ,M + Q) with pQ being 
the location of spinon excitation as 



M 



Q 

N 



N 

Do 
-Da 



(Toip)dp, 



Poiq)dq. 



(A-6) 
(A-7) 



Jn+l — Jn. — ^ + Sn 



(A-15) 



Equations (A-4)-(A-7]_lfiad to the ground-state distribu- 
tions (To{p) and po{q)^ 



<^o{p) = 



Poil) 




for \p\ < (1 — n)T^, 

for (1 - n)Tr < \p\ < (1 - fh)Tr, (A-. 

otherwise, 



where J^ £ [~{M + Q)/2, {M + Q)/2]. The position of 
a "hole" in the J^-configuration is denoted by J. When 
-(M + Q)/2 < J < /i or Iq+i < J < {M + Q)/2, 
a spinon is well-defined in the pseudomomentum space. 
The distribution functions cr(p) and p{q) must satisfy the 
following equations: 



aip)^e{B^\p\ 



1 

2^ 



D 



dqp{q)6{p - q) 



D 



l/(27r), 
0, 



for |g| < (1 — n)7r, 
otherwise. 



(A-9) 



We then obtain the ground-state energy per site 

-2 1 



fo 

t 



dpao{p)e{p) + —-{1-n) 



3n^ + m - 



-mh 



N 



dpa{p)6{p — p 



5{p~p^)Q{\p\~D) 



(A-16) 



1 



^2V" -- • ' ^-iff?)^-fhK (A-10) 

in the thermodynamic limit .EJ'cd) The ground-state mo- 
mentum is given by 



(A-11) 



qo ^ N dp(Jo{p){p + vr) 



p{q) = e{D - \q\) f dpa{p)5{q-p), (A-17) 

J -B 

where p^ is the value of p for p ~ /ig . Note that these 
equations fulfill the normalization conditions: 

(-B 

(T{p)dp 



!^(2 
2 ^ 



to). 



/ P{q)dq = 

J-D 



M + Q 


N 


Q + l 



N 



Now that we have obtained information on the ground 
state, the next task is to derive the dispersion relations 
of quasi-particles. We take the following steps: 

(i) We first give the configurations of J^ and li so that 
the quasi-particle can be created in the pseudomomen- 
tum space. 

(ii) We find distribution functions (j(j)) and p{q) under 
the configurations given in (i). 

(iii) In general, the distribution functions determining 
the excitation energy can be written in the form 

cTip)=ao{p) + il/N)ai{p), (A-12) 

p{q)=Po{q) + {l/N)p,{q). (A-IS) 



By introducing ai{p) by a{p) = ao{p) 
wc find 

1 



(A-lS) 

(A.19) 

(l/iV)ai(p), 



cri(p) = -^^ip-Ph)&i\p\ - (1 -n)7r) 
xe((l-m)7r- \p\). 



(A-20) 



We have used B = Bq and D = Dq in the large- iV limit. 
Then the excitation energy from the chemical potential 
is given by 



/ dp(Ji{p)e{p) - - 
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1.2 2^ h 



(A-21) 



as follows: 



where (1 — n)7r < |ph| < (1 — m)TT. For the associated 
momentum transfer — fc, we obtain 



dp<7i{p){p + tt) 



^{Ph+TT) 



(A-22) 



Therefore we obtain eq. (3.1). The constraint on ph re- 
sults in the following range of allowed momentum of up 
spinous: Tjm/2 < k < n7r/2 and (1 — n/2)7r < k < 
{l-m/2)TT. 

Next, in order to obtain the down-spinon dispersion, 
we consider the spin excitation in the subspace with 
{N, Q + 1, M), where an up spin is removed. The num- 
bers li (« = 1, 2, • • • , Q + 1) are chosen as 



/,; 



-Q/2 + i-l, 



(A-23) 



and the configuration of J^ (/x = 1, 2, • • • , Af + Q + 1) are 

^p+i-'^p = 1 + '5m,mo, (A-24) 

where J^ G [-(M + Q + l)/2, [M + Q + l)/2]. In a 
manner similar to that for the up spinon, we get the 
same expression of (Ti{p) as given by eq. (A-20). The 
second term of eq. (A-21) is replaced by +h/2 because 
the subspace we consider is different from that for the 
up spinon. 

A. 2 Antispinons 

In order to obtain the spectrum of antispinons, we 
consider the excitation in the subspace with {N, Q,M + 
1), where an up-spin flips into a down-spin. The numbers 
/j (i = 1, 2, • ■ • , Q) and J^ {p, = 1,2,- ■ ■ ,M + Q + 1) are 
chosen as follows: 

/, ==-Q/2 + z-l, (A-25) 

J^ = ~{M + Q)/2 + p, 

for ^= 1,2,- ••,M + Q, (A-26) 

Jm+q+1 - Jm+q > 1, (A-27) 

or 

/, = -g/2 + i, (A-28) 

J^ = -(M + Q)/2 + ^-2, 

for ^ = 2, 3, • • • , M + g + 1, (A-29) 

J2 - Ji > 1. (A-30) 

In the former (latter) case, p-p is defined as the value of 
p for J^ ~ Jm+q+1 (Ji)- The distribution function a(j>) 
can be written in the form 

o'(p) = o-o(rt©((l - ™)7r - IpI) 
1 



N 



Sip - pp)e{\p\ - [1 ~ m)Tr) 



xQ(tt-\p\), 



(A-31) 



Note that this equation fulfills the normalization condi- 
tion: / (j{p)dp = (M + Q + l)/N. We can identify ai{p) 



c^i(p) ^ 5{P - Pp)^{\p\ - (1 -m)iT) 



xe(7r- 


P\)- 


Then the excitation energy 


is given by 


^ = \ivl 


-n^) + h, 



(A-32) 



(A-33) 

and the associated momentum transfer is fc = pp + tt. 
Therefore we obtain the antispinon dispersion given by 
eq. (3.3). 

A. 3 Holons 

We consider the charge excitation in the subspace with 
{N, Q + 1, M — 1), where a down-spin is removed from 
the ground state. The numbers J^ (a* = 1, 2, ■ • • , M + Q) 
are chosen as follows: 



Ju 



-(M + Q)/2 + /i-l, 



(A-34) 



which means that the spin excitation (up-spinon exci- 
tation) takes the minimum energy that vanishes in the 
thermodynamic limit. On the other hand we choose li 

(i = 1, 2, • • • , Q -(- 1) either as 

/, = -{Q - l)/2 + i-l, for i = 1, 2, • • • , Q, (A-35) 

/q+1 - /q > 1, (A-36) 

or as 

h^-{Q-l)/2 + i-2, 

for i = 2,3, •••,0 + 1, (A^37) 

h~h-> 1. (A^38) 

In the former case we denote q-p = 2'kIq^i/N, and in the 
latter case ^p = 2ttIi/N. The distribution function p{q) 
can be written in the form 

p{q) = ^Oiil-n)n-\q\) 

+ ^S{q-qp)e{\q\-il-n)n). (A^39) 

On the other hand, the distribution function (j{p) satis- 
fies the following equation: 



a{p) = e{B-\p\ 



1 /-^ 

— + / dqp{q)6{p - q) 



dp'a{p')5{p-p') 



(A^40) 



Note that these equations fulfill the normalization condi- 
tions: J p{q)dq = {Q + 1) /N and J a{p)dp^ {M + Q)/N. 
By introducing ai{p) by a{p) = cro(p) + (1/^)cti(p), 
we find 

ai(p) = i5(p-qp)e(H-(l-n)7r) 

xe((l-m)7r- IpI). (A^41) 

We have used B = Bq in the large- iV limit. Then the 
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excitation energy is given by 
1 



M 



2 ^2n 
TT 



(A-42) 



4V-P ' ■ y- 

The second term conies from tt^Q{\ — l/A^^)/3 of eq. 
([A-1| ). The momentum transfer is given by fc = {q^ + 
7r)/2. The relation between e and k leads to eq. (3.4). 
The momentum k is allowed for fhix jl < k < nix jl and 
(1 — n/2)n < k < (1 — m/2)7r. The lowest excitation 
energy, which occurs at fc = n7r/2 or (1 — n/2)7r, cor- 
respforuis to the absolute value of the chemical poten- 

tiaiee 

A.J^ Antiholons 

In order to obtain the antiholon dispersion, we con- 
sider the excitation where two electrons are added. We 
assume that one of the two electrons has up-spin and 
the other has down-spin. Then the subspace becomes 
{N,Q-2,M+1). The numbers of J^ {n = l,2,---,M-h 
Q — 1) are chosen as 



Jp = -(M + Q)/2 + M, 
and the configurations of /^ (i = 1, 2, • • • 

li+l — li = ^ + ^i,kn 



(A-43) 



2) are 



(A-44) 



where /, e [-{Q - 2)/2, (Q - 2)/2]. The posi_tion of a 
"hole" in the /^-configuration is denoted by /. When 
-{Q - 2)/2 < / < (Q - 2)/2, an antiholon is well- 
defined in the pseudomomentum space. The distribution 
function p{q) can be written in the form 



Pi'l) = 



i l^('^ - '^'^^ 



e((l-n)7r-|g|), (A-45) 



where q^ is the value of q for i = iq. Note that this 
equation fulfills the normalization condition: J p{q)dq = 
[Q - 2)/7V. 

By introducing (Ji{p) by (j{p) == cro(p) + {^/N)(Ti{p), 
we find 



ai(p)--<5(p-(Zh)e((l-nV-b|). 
Then the excitation energy is given by 



\i<^l 



r 



(A-46) 



(A-47) 



The second term comes from 7r^Q(l — l/A^^)/3 of eq. 
(A-1). The momentum transfer is given by —k = —{qh + 
tt). Therefore we obtain eq. (3.5). 
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Table I. List of all the possible excitations from the ground state 
with finite magnetization. 



Fig. 8. S'^^ {q,u!)/4 of the l/r^-type supcrsymmctric f- J model 
with 16 sites, 2 holes and (a) m = 0.25 and (b) m = 0.5. The 
broken line indicates the magnetic field h and its double 2h. The 
hatched area indicates the region where holons and antiholons 
do not contribute at all. 



Fig. 1. Dispersion relations of (a) up and down spinons, (b) anti- 
spinons, (c) holons and (d) antiholons as given by eqs. (3.1)-(3.5). 
Allowed ranges of momenta and the right and left branches are 
shown for each species. 



Fig. 9. N(q,u!) of the 1/r^-type supcrsymmetric t- J model with 
16 sites, 2 holes and (a) rh = 0.25 and (b) m = 0.5. The broken 
line indicates the magnetic field h. The hatched area indicates 
the region where spinons and antispinons do not contribute at 

all. 



Fig. 2. Skew Young diagram corresponding to the motif 
010|10|110|10|100. There are other possibilities to put 1,2 or 
o than the one shown in this figure. 



Fig. 10. S' 
16 sites, n 



^{q) of the l/r^-type supersymmetric i- J model with 
= 0.25 and two different values of filling. 



Fig. 3. S ^{q,Lu)/4 of the Haldane-Shastry model with 16 sites 
and (a) fh = 0.25 and (b) m = 0.5. The intensity of each pole is 
proportional to the area of the circle. The solid lines show the 
dispersion lines of quasi-particles in the thermodynamic limit. 
This way of representation applies to all figures to follow. 



Fig. 11. S^^{q) of the 1/r^-type supersymmetric 4- J model with 
16 sites, fh = 0.25 and two different values of filling. 



Fig. 12. ^{q) of the 1/r'^-type supersymmetric t-J model with 
16 sites, n = 0.875 and two different values of magnetization. 



Fig. 4. S ^{q,Lo)/4 of the 1/r^-type supersymmetric i- J model 
with 16 sites, 2 holes and (a) fh = 0.25 and (b) fh = 0.5. The 
hatched area indicates the region where only two spinons con- 
tribute. 



Fig. 13. S^^{q) of the nearest-neighbor-type t-J model with 16 
sites, fh = 0.25 and two different values of filling, (a) J/t = 2; 
(b) J/t = 0.5. 



Fig. 5. S^^{q, uj) of the Haldane-Shastry model with 16 sites and 
(a) fh = 0.25 and (b) fh = 0.5. The broken line indicates the 
magnetic field h. 



Fig. 14. 5'^^((jr) of the nearost-neighbor-type t-J model with 16 
sites, fh = 0.25 and two different values of filling, (a) J/t = 2; 
(b) J/t = 0.5. 



Fig. 6. 5^^ (g, (i)) of the 1/r^-type supersymmetric i- J model with Fig. 15. 

16 sites, 2 holes and (a) fh = 0.25 and (b) fh = 0.5. The broken sites, 

line indicates the magnetic field h. The hatched area indicates J/t = 
the region where holons and antiholons do not contribute at all. 



N{q) of the nearest-neighbor-type t-J model with 16 
n = 0.875 and two different values of magnetization, (a) 
2; (b) J/t = 0.5. 



Fig. 7. S'^^ {q,Lo)/4 of the Haldane-Shastry model with 16 sites 
and (a) fh = 0.25 and (b) fh = 0.5. The broken line indicates 
the magnetic field h and its double 2h. 



This figure "figurel.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/0009075vl 



This figure "figurell.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/0009075vl 



Table 1 : List of all the possible excitations from the ground state with finite 
magnetization. 



Dynamical quantities 


Excitation contents 


S-+{q,u) 


(sl, hL) + h + (hR, sr) 




2sr, (L ^ R) 


S^^{q,uj) 


(sl, hL) + h + (hR, sr) + sr, (L ^ R) 




(sL, hL) + h + (hR, Sr) 




2Sr + Sr, (L ^ R) 




2Sl + Sr, (L ^ R) 




2sr, (L ^ R) 


S+-{q,u) 


Sl + (sl, hL) + h + (hR, sr) + sr 




(sl, hL) + h + (hR, sr) + SR, (L ^ R) 




(sl, hL) + h + (hR, Sr) 




Sl + 2sR + SR, (L ^ R) 




2sr + Sr, (L ^ R) 




2Sl + Sr, (L ^ R) 




2sr, (L ^ R) 




Sl, (L ^ R) 


N{q,iu) 


(sl, hL) + h + (hR, Sr) + Sr, (L ^ R) 




(sl, hL) + h + (hR, Sr) 




h + 2hR, (L ^ R) 



